Neutron star tidal disruption in mixed binaries: 
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We study the tidal disruption of neutron stars in black hole-neutron star coalescing binaries. We 
calculate the critical orbital separation at which the star is disrupted by the black hole tidal field 
for several equations of state describing the matter inside the neutron star, and for a large set of the 
binary parameters. When the disruption occurs before the star reaches the innermost stable circular 
orbit, the gravitational wave signal emitted by the system is expected to exhibit a cutoff frequency 
vcwtide, which is a distinctive feature of the waveform. We evaluate vawtide and show that, if this 
frequency will be found in a detected gravitational wave, it will allow one to determine the neutron 
star radius with an error of a few percent, providing valuable information on the behaviour of matter 
in the stellar core. 
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I. INTRODUCTION 

The coalescence of neutron star-neutron star (NS-NS) 
and black hole-neutron star (BH-NS) binaries is one of 
the most promising sources of gravitational waves to 
be detected by ground based gravitational wave detec- 
tors like Virgo and LIGO [l|. These detectors have 
now reached the planned sensitivity and they will evolve 
toward a second generation, the advanced (Virgo and 
LIGO) detectors, with a sensitivity enhanced by an or- 
der of magnitude. Furthermore, a design study for an 
even more sensitive third generation detector, ET (Ein- 
stein Telescope), is in progress @. 

In a recent study, based on a population synthesis ap- 
proach the formation and the evolution of compact 
binary systems has been followed from the onset of star 
formation, both in the Galactic field, where their massive 
binary progenitors evolve in isolation, and in dense clus- 
ters, where they can form at high rates due to dynamical 
interactions. The authors estimate that advanced LIGO 
should detect the merger of NS-NS binaries at a rate of 
~ 15 events per year and the merger of BH-NS binaries 
at a rate of ~ 1 per year. Similar estimates hold for the 
advanced version of Virgo. Thus, it is reasonable to ex- 
pect that in a near future we may be able to detect the 
gravitational wave (GW) signals emitted by these sources 
and study their features. 

A further reason to be interested in these coalescing bi- 
naries is that they have been proposed as providing the 
engine for short gamma-ray bursts (SGRBs); the detec- 
tion of a gravitational wave signal emitted by one of these 
systems in coincidence with a SGRB would validate this 
model, thus clarifying one of the most interesting open 
issues in astrophysics. 

In this paper we study the disruption of a neutron star 
in a BH-NS coalescing binary, to envisage a method to 
extract information on the equation of state (EOS) of 
matter in the NS interior from the detection of the grav- 
itational wave signal emitted in the process. We use the 



affine model approach 0)0,05 which treats the NS 
as an extended object, which responds to its self-gravity, 
to its internal pressure forces and to the relativistic BH 
tidal field, under the constraint that its shape is always 
that of an ellipsoid. In the original formulation of this 
model, the internal structure of the star was treated at 
a Newtonian level and the EOS was assumed to be poly- 
tropic. In [8] (to be referred to as Paper I hereafter) 
we have improved this approach by introducing a gen- 
eral relativistic description of the stellar structure, and 
generalizing the relevant equations to include more gen- 
eral equations of state. This improved approach has been 
applied to study quasi-equilibrium sequences of BH-NS 
binaries, and to determine the critical orbital separation 
at which the star is torn apart by the black hole tidal 
field. 

When the NS is disrupted before reaching the Inner- 
most Stable Circular Orbit (ISCO), the emitted GW sig- 
nal is expected to change abruptly and its amplitude is 
expected to decrease sharply; such signal should exhibit 
a cutoff frequency vowude, corresponding to the orbit at 
which the disruption occurs, i.e. VGWtide — ^v rb tide- 
This prediction is confirmed by numerical simulations 
where the coalescence of a black hole and a neutron star 
and the tidal disruption of the star have been studied i) 
in the framework of Newtonian gravity @, [T(J, and ii) 
in full general relativity [ll|, fl2| . 

The frequency cutoff is a distinctive feature of the 
waveform emitted by a BH-NS coalescing binary and in- 
dicates the disruption of the star. This was first pointed 
out in [l3j , where a relation between vcwude and the stel- 
lar radius was derived, using the formulation of the affine 
model appearing in [14| and describing the inspiralling 
of an incompressible, homogeneous, Newtonian ellipsoid 
moving on circular orbits around a rotating black hole. 
Hints on the role of the compressibility were derived in 
using some results of [6] , where the same process was 
studied for a polytropic, Newtonian star. 

In this paper we study the relation between the cutoff 
frequency and the neutron star EOS. Using the same ap- 
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proach developed in Paper I, we explicitly compute for 
which values of the masses of the binary components and 
of the black hole angular momentum the star is tidally 
disrupted before reaching the ISCO, for a variety of real- 
istic equations of state. Knowing the disruption distance 
rude, we evaluate fGWtide and show how a measure of 
this quantity from a detected GW signal may be used to 
infer interesting information on the equation of state of 
matter in the neutron star interior. 

It is known that the fully relativistic numerical study 
of the last phases of BH-NS binary coalescence is a quite 
difficult task, especially for large values of the mass ra- 
tio q = M BH /M N s', this is due to a number of rea- 
sons, which include the lack of symmetry, the difficulty in 
evolving the NS while handling the BH singularity, and 
the prohibitively high computational costs required to 
span the parameter space (mass ratio, EOS, spins etc.). 
For these reasons, the literature on the subject is limited 
to a restricted number of studies, most of which are in 
the framework of Newtonian gravity. BH-NS coalescence 
in general relativity has been studied in the following 
papers (see [l5[ for a recent review on the subject): non 
spinning black hole and Mbh » Mns in H6]; q — 10 in 
q = 1,3,5 infill; and M BH = i _3.2M , M BH = 4M Q 



and M NS = 1.4M Q in QJ]. In [T^, quasi-equilibrium 
sequences of black hole-neutron star binaries have been 
studied in general relativity for q = 1, 2, 3, 5, 10. The case 
of a rotating black hole has been considered in ll| , where 
the coalescence has been studied for three values of the 
black hole angular momentum (cl/Mbh = —0.5,0,0.75) 
and for q = 3, whereas in [l2| the case q = 10 has been in- 
vestigated. It should be stressed that in all these studies 
the neutron star is modeled using a polytropic equation 
of state. 

The improved affine model approach which we use in 
our study is an approximate method, but it has the ad- 
vantage of allowing one to explore a large region of the 
parameter space, including large values of the mass ratio 
q = Mbh/Mns, an d to use modern equations of state 
to model matter in the neutron star interior, at a much 
lower computational cost than fully relativistic simula- 
tions. 

With respect to ref. [lU our study introduces several 
novelties: we describe the internal structure of the star 
using general relativity, we explicitly take into account 
the stellar compressibility and describe the matter inside 
the star using realistic equations of state; in addition, in 
[3 vcwtide was provided by the Kerr geometry in the 
approximation of a point mass NS following geodesies in 
this geometry, while we determine this quantity using the 
2.5 Post Newtonian equations of motion for the orbital 
dynamics of a binary system. 

The plan of the paper is the following. In Section |TT] 
we briefly describe the affine approach; in Section 11111 we 
discuss how to determine the cutoff frequency vawtide- 
The equations of state we employ are briefly described 
in Section HVt in Section [V] we discuss the errors which 
affect the measure of the binary parameters from a de- 



tected GW signal, and how they propagate and affect 
the evaluation of vowude and of the neutron star radius. 
The results of our study are reported in Section IVII and 
in Section [VIII we draw our conclusions. 



II. THE IMPROVED AFFINE MODEL 
EQUATIONS 

The equations of the improved affine model approach, 
which we need in order to determine the radial distance 
Tude at which the neutron star is torn apart by the tidal 
interaction with the black hole, are described in detail in 
Paper I. Here we will summarise the assumptions under- 
lying this approach and we shall write only the equations 
to be solved. 



1. Relevant assumptions and equations 

We consider a star in the tidal field of a Kerr BH 
whose center of mass follows equatorial, circular orbits; 
while moving, the star maintains an ellipsoidal shape; 
more precisely it is a Riemann-S type ellipsoid, i.e. its 
spin and vorticity are parallel and their ratio is constant 
(see [13|). The NS equilibrium structure is determined 
using the stellar structure equations of General Relativ- 
ity, while dynamical behaviour is governed by Newtonian 
hydrodynamics improved by the use of an effective rela- 
tivistic self-gravity potential. The equations for the NS 
deformations are written in the principal frame, i.e. the 
frame associated with the principal axes of the stellar el- 
lipsoid. Tidal effects on the orbital motion are neglected, 
as well as the perturbation that the star induces on the 
BH. 

We study the evolution of the system in the quasi- 
equilibrium approximation, i.e. we neglect all time 
derivatives appearing in the ordinary differential equa- 
tions of the model. Physically this means that we assume 
that the NS follows a quasi-equilibrium sequence during 
the BH-NS coalescence; according to the model, more- 
over, the circulation C of the fluid along this sequence is 
constant. In particular, we set C = 0, that is, we con- 
sider the NS fluid to be irrotational. In the principal 
frame, the fluid variables of the improved affine model 
are five: the three principal axes of the stellar ellipsoid 
Oi, &2 and 03, the angular frequency of the internal fluid 
motion A and the star spin Q. The axes a\ and a 2 both 
belong to the orbital plane, while 0,3 is perpendicular to 
it; di indicates the axis that lies along the binary orbital 
radius. The quasi-equilibrium equations governing these 
variables are: 
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where V and .M are, respectively, the effective relativistic 
self-gravity potential and the scalar quadrupole moment 
of the NS at spherical equilibrium (see below). Rns is 
the NS radius and the A-'s are defined as 
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The Cjj's denote the components of the BH tidal tensor 
in the principal frame; <j> — which is connected to fl by 
4> = il — is the angle that brings the parallel-propagated 
frame in the principal frame, by a rotation around the 03 
axis [2lJ . Finally, \& is an angle that governs the rotation 
of the parallel-propagated frame in order to preserve its 
parallel transport, and its time evolution is given by 
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where Mbh an d r are the BH mass and the BH-NS or- 
bital separation respectively. Notice that the last equa- 
tion and Eq. (j4j imply that 
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The fifth equation for the five fluid variables is provided 
by the definition of the circulation 
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so that for irrotational fluids one has 
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The effective relativistic self-gravity potential is given 
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where d&xov /df is given by the Tolman-Oppenheimer- 
Volkoff (TOV) stellar structure equations 
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The scalar quadrupole moment A4 is defined as 
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and, like V, must be calculated at spherical equilibrium. 
Finally, the relevant Kerr BH tidal tensor components 
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a is the black hole spin parameter, and E and L z are, re- 
spectively, the energy and the z-orbital angular momen- 
tum per unit mass. Since we consider circular equatorial 
geodesies, E and L z are 
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2. Numerical integration 

We solve Eqs. with the aid of Eqs. Q, © 

and of definitions (|5]l, ([TO]). ([TT]), (fT2" ]) -([T3 1) . by adopting 
a multidimensional Newton-Raphson scheme (22l | in or- 
der to determine the values of the axes of the ellipsoid 
for each quasi-stationary orbit, identified by the orbital 
separation r. 

We start by solving the TOV stellar structure equa- 
tions for a non-rotating spherical neutron star in equi- 
librium. We then fix the black hole spin parameter a 
and the binary mass ratio q and place the star at a dis- 
tance ro ^> Rns from the black hole (we obviously make 
sure that the sequence we obtain is independent of ro). 
Subsequently, we gradually reduce the orbital separation 
and solve Eqs. ([I])-© at each step and monitor the star 
axes until a critical separation rude is reached, at which 
a quasi-equilibrium configuration is no longer possible. 
This critical distance physically corresponds to the tidal 
disruption of the neutron star. It may be identified by 
exploiting the fact that the Newton-Raphson algorithm 
cannot find any solution to the system of equations, or by 
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calculating numerically the derivative dr norm / 'd(a,2/ '«i), 
where 

norm ~ RnsKMbh) ' 

and keeping track of it since it tends to zero at tidal dis- 
ruption: both methods yield the same values of rude- We 
mention that an alternative approach to evaluate rude ln 
a Newtonian framework, based on the estimate of the 
Roche lobe radius, has been used in (23|. If the tidal 
disruption is not encountered, the quasi-equilibrium ends 
when the neutron star surface crosses the black hole hori- 
zon and hence the coalescence terminates with a plunge. 

The quantity rude has to be compared with the value 
of the radius of the innermost stable circular orbit risen, 
which is determined by using the formulae derived in 24] 
for a point mass in the gravitational field of a Kerr BH: 
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where the upper (lower) sign holds for co- (counter- 
rotating orbits. We remind the reader that if rude > 
fisco the star is disrupted before the merger starts, and 
the gravitational signal will exhibit a cutoff at a frequency 

VGW tide- 



Ill. DETERMINATION OF THE CUTOFF 
FREQUENCY vowude 

To compute VGWtide, we model the inspiral of the 
mixed binary by means of a post-Newtonian (PN) ap- 
proach and truncate the inspiral when the orbital sepa- 
ration reaches rude] we then read off the orbital frequency 
at the tidal disruption v or b tide- this is related to the GW 
frequency by v G wude = 2v orb tide- 

We follow a Hamiltonian approach. The conservative 
part of the two-body Hamiltonian is known up to order 
3PN, e.g. [25|; however, we shall truncate it at 2PN order 
since we will use GW dissipation terms of order 2.5PN 
[26| . The PN-expanded Hamiltonian for the relative mo- 
tion of the BH-NS binary is 
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where Gn is the gravitational constant, c is the speed 
of light, Mrot — M B h + Mns is the total mass of 
the system, /i = M BB M^slMTot is its reduced mass, 
rj = n/MTot is the symmetric mass ratio, P r and P v are 
the conjugate variables of the orbital separation r and the 
orbit angle coordinate ip, L = r x P is the orbital angular 
momentum and Jbh is the spin angular momentum of 
the black hole. The vectors L and Jbh are both perpen- 
dicular to the orbit plane; the magnitude of the latter is 
equal to clMbh and such vector is not evolved as a dy- 
namic variable in order to follow the spirit of the afiine 
model, according to which the presence of the NS does 
not influence the BH. The spin-orbit (SO) and spin-spin 
(SS) contributions to therefore reduce to 
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The dynamics described by this Hamiltonian is sup- 
plied with GW dissipation by means of the 2.5PN non- 
conservative terms 
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to be included in the Hamilton equations of motion as 

fr (31) 

U (32) 
F r (33) 
■F v . (34) 

Within this Hamiltonian approach, once the binary or- 
bit is evolved until r = rude, the GW cutoff frequency 
due to the NS tidal disruption vowude is given by 
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IV. EQUATIONS OF STATE 

The equation of state of matter at densities larger than 
the saturation density of nuclear matter, po — 2.67- 10 14 
g/cm 3 , is uncertain. At these densities, which are typical 
of a neutron star core, neutrons can no longer be con- 
sidered as non interacting particles, and different ways 
of modeling their interactions lead to a different com- 
position, which may include heavy baryons, quarks etc.. 
All available EOSs of strongly interacting matter have 
been obtained within models based on the theoretical 
knowledge of the underlying dynamics and constrained, 
as much as possible, by empirical data. Modern EOSs 
are derived within two main, different approaches: non- 
rclativistic nuclear many-body theory (NMBT) and rel- 
ativistic mean field theory (RMFT). As representative 
of NMBT we choose two EOSs named APR2 and BBS1. 
For APR2 matter consists of neutrons, protons, electrons 
and muons (n,p, e, p) in weak equilibrium. The Hamilto- 
nian includes two- and three- nucleon interaction terms; 
the two-nucleon term is the Argonne t>x8 potential 27J, 
the three- nucleon term is the Urbana IX potential 28|. 
The many-body Schrodinger equation is solved using a 
variational approach (29l |30| . The calculations include 
rclativistic corrections to the two-nucleon potential, aris- 
ing from the boost to a frame in which the total mo- 
mentum of the interacting pair is non- vanishing. These 
corrections are necessary to use phenomenological poten- 
tials, describing interactions between nucleons in their 
center of mass frame, in a locally inertial frame associ- 
ated with the star. The maximum mass for this EOS is 
M max = 2.2OM . 

The matter composition of BBS1 is the same as in 
the APR2 model. The EOS is obtained using a slightly 



different Hamiltonian, including the Argonne vig two- 
nucleon potential and the Urbana VII three-nucleon po- 
tential [31]. The ground state energy is calculated using 
G-matrix perturbation theory [32j . The maximum mass 
is M max = 2.OlAf . 

As representative of the RMFT approach, we choose 
the EOSs named BGN1H1, GNH3 and BPAL. 

The Balberg-Gal (BGN1H1) EOS 33] describes mat- 
ter consisting of neutrons, protons, electrons, muons and 
hyperons (£, A and S) in equilibrium. The effective po- 
tential parameters are tuned in order to reproduce the 
properties of nuclei and hypernuclei according to high 
energy experiments. This EOS is a generalization of the 
Lattimer-Swesty EOS [34|, which does not include hy- 
perons. The maximum mass is M max = 1.63 Mq. 

The Glendenning (GNH3) EOS [HI considers n,p,e,/j, 
up to a certain density pn — 2po! beyond this point, 
additional baryon states (such as the A and the hyperons 
A, S, S) and the mesons n, a, p, uj, K, K* are introduced. 
Below the hadronization density pn the EOS is very stiff 
but causal; the appearance of hyperons strongly softens 
the EOS because they are more massive than nucleons 
and when they start to fill their Fermi sea they are slow 
and replace the highest energy nucleons. The maximum 
mass is M max = 1.96 Mq. 

The three EOSs BPAL12, BPAL22 and BPAL32, are 
derived using a density dependent nucleon-nucleon ef- 
fective interaction (as for Skyrme nuclear interactions) 
[3f|. Matter is composed of n,p,e,p in weak equilib- 
rium, no hyperons are present and the EOS parameters 
are fixed to reproduce the saturation properties of nu- 
clear matter. BPAL12 is particularly soft, with a nu- 
clear incompressibility k = 120 MeV. Since the empirical 
value commonly accepted ranges within k 6 (220 — 270) 
MeV, it is clear that BPAL12 has to be considered as 
an EOS soft extreme, but still compatible with astro- 
physical observations, since the maximum mass it pre- 
dicts is M max = 1.45 M . BPAL22 and BPAL32 are 
two different versions of the BPAL EOS, correspond- 
ing to more realistic values of k: k — 180 MeV for 
BPAL22, with a maximum mass M max = 1.72 Mq, 
and k = 240 MeV for BPAL32, with a maximum mass 
M max = 1.93 Mq. At a density of about one half of the 
nuclear saturation density we match the EOS of the core 
with a crust, which is composed of three layers: the inner 
layer is the Douchin-Haensel (SLy4) crust [33, [38|; for 
10 8 < p < 10 11 g/cm 3 and for p < 10 8 g/cm 3 , we use, re- 
spectively the Haensel-Pichon (HP94) EOS [39| . and the 
Baym-Pethick-Sutherland (BPS) EOS 03 . 

In Fig. [T] we show the mass-radius diagram for the 
selected EOSs. They clearly exhibit very different be- 
haviours. To some extent, BPAL12 and GNH3 can 
be considered, respectively, as soft and hard extremes. 
APR2 has a wide mass interval where the radius is al- 
most insensitive to mass variations; BGN1H1 also shows 
this feature, albeit for a more restricted mass interval, 
whereas the BPAL and GNH3 EOSs do not. BGN1H1, 
on the other hand, has a sudden softening which does 
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FIG. 1: ( Color online) The neutron star mass is plotted versus 
the stellar radius for the EOSs used in this paper. 



not characterise the other three equations of state; this 
softening is due to the appearance of hyperons in the core 
of the neutron star above a critical central density. 



V. ERROR BAR ON R NS AND ON vawtide 

Let us suppose that in a BH-NS coalescence the neu- 
tron star is disrupted before reaching the ISCO, and that 
the gravitational wave interferometers detect the emit- 
ted signal which, as discussed in previous sections, has 
the form of a chirp truncated at the frequency VGWMde- 
With a suitable data analysis, the values of the symmet- 
ric mass ratio rj = Mbh^ns/ (Mbh + M^s) 2 and of the 
chirp mass M. = t) 3 ' 5 (Mbh + Mns) can be extracted 
from the data with a certain error. These errors have 
been evaluated for non-spinning c omp act binary sources, 
including up to 3.5PN terms, in |41| . [42{ ■ For instance, 
they are displayed in Fig. 12 of [42|], which shows the 
one-sigma fractional errors in A4 and r\ for non-spinning 
binary black hole sources as a function of the total mass 
of the system. These errors reduce when the dynamical 
evolution of spins is included, so that they may be re- 
garded as upper limits. Since a similar analysis has not 
been performed in the BH-NS case, in what follows we 
shall adopt these errors as appropriate also for mixed bi- 
naries. In the following we shall consider neutron stars 
whose mass is in the range 1.2 — 1.6 solar masses and 
mass ratios from 3 to 15, therefore the total mass will 
range within ^5 — 26 M©. In this region the fractional 
error on Ai is smaller than ~ 10~ 3 and of the order of 
~ 1 — 3 ■ 10~ 2 for rj. These data refer to advanced LIGO 
assuming the source at a fixed distance of 300 Mpc. 

Since in our case the mass parameters of the binary 
are the neutron star mass and the mass ratio, we express 



TABLE I: Absolute errors on the neutron star mass and on 
the mass ratio. Aq is independent of the neutron star mass. 
AMns is expressed in solar mass units. 
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By propagating the errors we find the following expres- 
sions for the one-sigma fractional errors in the neutron 
star mass and in the mass ratio: 
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The absolute error on M^s and q for the binaries we 
consider are given in Table HI We shall now discuss how 
the errors which affect the estimate of the binary param- 
eters influence the evaluation of the neutron star radius 
and of vawtide- Since we do not know the error on the 
black hole angular momentum, in what follows we shall 
assume that, with an accurate post-detection data analy- 
sis, a/MsH could be measured with a 10% accuracy. As 
an example, let us suppose that the neutron star mass 
and the mass ratio measured from a detected signal are, 
say, 1.4 Mq and 5, respectively, plus or minus the corre- 
sponding error which can be found in Table [I] In Fig. [2] 
we plot the neutron star radius versus VGWtide, evaluated 
as explained in section UTTl for any possible combination 
of the following data: 



M NS = (1.4±0.02)M ( 



©> 



M 



BH 



0.5±0.05, q = 5±0.15, 
(39) 



and for a given EOS, for instance BPAL12. 



The data are spread in a region which can be identified 
by a vertical and a horizontal error bar, indicated as a red 
cross in the figure: the vertical error corresponds to the 
maximum and minimum values of the neutron star mass. 
The horizontal error is identified by the points (Mns — 
AM N s,q+Aq,a min ) and (M N s + AM NS ,q~ &q,a ma x)- 
In the following figures, where we will show the graphs 
of Rns versus VGWtide for different values of the binary 
parameters and for different EOSs, we shall plot directly 
the cross, determined as in Fig. [5J 
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EOS BPAL12 



q+Aq, a m 



M NS- AM NS 
q-Aq, a max 



M NS +AM NS 
q+Aq, a min 



M NS +AM NS 
q-Aq, a max 



1400 1450 
v GWtide ( Hz ) 



FIG. 2: (Color online) The NS radius is plotted versus 
VGWtide for any possible combination of the following data: 
Mns = 1.4 ± 0.02 Mq, a/MsH = 0.5 ± 0.05, q = 5 ± 0.15 , 
and for the EOS BPAL12. The horizontal and vertical sp 
of the data identify the error bars (red cross). 



VI. RESULTS 

The results of our calculations are summarised in the 
Figs. HH 



In Fig. |3]we plot the NS radius, Rns, versus the cutoff 
frequency, VGWtide, for a black hole with angular mo- 
mentum cl/Mbh — 0.9 ± 0.09, and for a neutron star 
with mass M NS = 1.4 Mq ± AM NS in Fig. E^a), and 
M NS = 1-2 M Q ± AM NS in Fig. EJb). The data are 
plotted for different values of the mass ratio q ± Aq (we 
omit writing "±Aq" in the figures) and for the consid- 
ered EOSs. The values of AAIns and Aq are given in 
Table HI For each value of q, the continuous black lines 
are parabolic fits of the data corresponding to each EOS. 
Fig.[3]shows that, for a given EOS, as q increases VGWtide 
decreases. Moreover, for q = 15 the data corresponding 
to the EOSs BBS1, BPAL22, APR2, BPAL12 are miss- 
ing. This behaviour is easily understood if we plot the 
radius at which disruption occurs, r t ^ e , versus VGWtide 
for different values of q. For instance, in Fig. 0] this plot 
is done for M^s — 1-4 Mq and ajMsH = 0.9 for the 
EOSs BGN1H1 and BPAL12. For comparison, we also 
plot the value of rjsco, given by Eqs. (ITS)) , versus the 
frequency vqw isco of the gravitational signal emitted 
when a point mass of mass M = M^s, orbiting around 
a Kerr black hole of mass Mbh = qM and angular mo- 
mentum cl/Mbh = 0.9, reaches rrsco- For the BGN1H1 
star (upper curve) we see that r t ide is larger than rjsco 
for all values of q; thus the GW signal emitted by these 
systems will exhibit a frequency cutoff at VGWtide- More- 
over, as q increases the value of r t ide increases, i.e. the 
star is disrupted at larger distances from the black hole. 
As a consequence vcwude is a decreasing function of q, 



q=15 q=10 




a)M NS =1.4 M sun ±AM NS 



1200 1400 
v GWtide ( Hz > 



q=15 q=10 




600 800 1000 1200 1400 1600 1800 20! 

v GWtide ( Hz ) 

FIG. 3: (Color online) Rns is plotted versus VGWtide, for a 
black hole with angular momentum a/ Mbh = 0.9 ± 0.09 and 
a neutron star with mass Mns = 1.4 Mq ± AMms (upper 
panel) and Mns = 1-2 Mq ±AMns (lower panel) . The data 
are plotted for different values of the mass ratio q ± Aq and 
for the EOSs considered in section UV\ The values of AMns 
and Aq are given in Table Black continuous lines are the 
parabolic fits of the data corresponding to each EOS at fixed 



as shown in Fig. [3] for all EOSs. 



A similar behaviour is shown by the EOS BPAL12 (mid- 
dle curve) up to q = 10. At that point the curve r t ide- 
versus-i'Gwtide crosses the lower curve r/sco-versus- 
v gw iSGOi and for larger q r t ide becomes smaller than 
rjsco', this means that, for q > 10, the BPAL12 star 
would be disrupted after reaching the ISCO and the emit- 
ted signal would not exhibit a frequency cutoff. For this 
reason the part of curve between the points correspond- 
ing to q = 10 and q = 15 is indicated as a dashed line. 
Going back to Fig. [3l the above discussion clarifies why 
for q = 15 the data for the EOSs BBS1, BPAL22, APR2 
and BPAL12 are missing: the star merges with the black 
hole without being disrupted, and there is no frequency 
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BGN1 H1 : r tide (v Gwtid8 ) « 








BPAL12:r tide (v Gwlide ) b " 
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ioVvl_____ 
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M NS= 14 M sun 









500 1000 1500 2000 2500 3000 3500 

v GWtide ( Hz ) 

FIG. 4: (Color online) The tidal disruption radius, rude, is 
plotted versus VGWtide for the EOSs BGN1H1 (upper curve) 
and BPAL12 (middle curve), assuming Mjsrs = 1.4 Mq and 
cl/Mbh = 0.9. The numbers label points corresponding to dif- 
ferent values of q. For the same values of q we also plot the 
value of risco (lower curve) versus the corresponding fre- 
quency vgw isco (see text). The dashed part of the curve 
for BPAL12 refers to values of q for which the star is not 
disrupted before the ISCO, i.e. for which rude < risco- 



a/M RH =0.9±10% 




600 800 1000 1200 1400 1600 1800 2000 



v GWtide ( Hz ) 

FIG. 5: (Color online) The same plot as in Fig. but for 
M N s = 1.6M Q ± AM NS . 



Fig. [T] moreover, the vertical error bars are larger, be- 
cause for this EOS Mns — 1.6M is very close to the 
maximum mass, where the mass-radius curve is nearly 
flat (see Fig. [J). 



cutoff in the emitted GW signal. 

In Fig. [3] we also notice that the vertical errors are al- 
ways much smaller than the horizontal ones, except for 
the EOS BPAL12, when M NS = 1.4M ± AM NS . The 
reason is that the vertical error is due to the error on 
the neutron star mass, which we recall is of the order of 
a few percent; the curves corresponding to the different 
EOSs in the mass-radius diagram of Fig.Q]show that, for 
M NS = 1.4 M Q or M NS = 1.2 M , the stellar radius 
does not change significantly for such a small change in 
Mns- However, Mns — 1-4 M Q is close to the maxi- 
mum mass of the BPAL12 EOS, and in its neighborhood 
the mass-radius curve is almost flat. For this EOS even 
a small change in the mass corresponds to a significant 
change in the radius. 

Fig. GJb) shows that, when M NS = 1.2M Q ± 
AM NS , for the pairs of EOSs (BBS1,BPAL22) and 
(APR2,BPAL12) the data nearly coincide. The reason 
is that, as shown in Fig. [TJ for Mns = 1-2 Mq the stars 
corresponding to these pairs of EOSs have nearly the 
same radius. For the same reason the values of vgw tide 
for BGN1H1 and BPAL32 are very close. 



In Fig. [5] the same quantities of Fig. [3] are plotted for a 
larger mass, Mns — 1-6M ± AMns- The data for the 
EOS BPAL12 are missing because the maximum mass 
allowed by this EOS is M max = 1.45 Mq. A compari- 
son with Fig. [3]Ja) shows that, for Mns = 1-6M , the 
data for the EOS BGN1H1 move at the bottom of the 
figure, i.e. for these masses the NS radius is smaller than 
that given by the other EOSs, as can also be seen from 



In Fig. [6] the same quantities of Fig. [3] are plotted 
for a black hole with angular momentum cl/Mbh 
O.o ± 0.05. Fig. H^a) refers to a neutron star with 
mass Mns = 1-4 Mq ± AMns, whereas in Fig. \&[b) 
M NS = 1-2 M ± AM N s- The data for q = 10 and 15 
are missing because no disruption occurs for the consid- 
ered EOS: tidal disruption is favoured by large values of 
the black hole angular momentum. 

If we compare this figure with Fig. [3] we see that, for 
a given mass and EOS, the value of the cutoff frequency 
is smaller if the black hole angular momentum decreases. 
This means that the tidal disruption radius rude increases 
as a decreases. It should be noted that for a fixed value 
of Mns and of q risco also increases as a decreases; it 
increases faster than rude, therefore there exists a critical 
value of the black hole angular momentum below which 
the star is not disrupted before the ISCO. In a similar 
way, for fixed a, Mns and EOS, there exists a critical 
value of q above which no disruption occurs. As an ex- 
ample, in Table [TT] we give the values of q ma x and the 
corresponding vowtide, which is the minimum value of 
the cutoff frequency to be expected, for the binaries con- 
sidered in Figs. [3] and [B] 



A direct comparison of the effect of the angular momen- 
tum on VGWtide is displayed in Fig. [7] where the plot of 
Rns versus v GWti de is done for M N s = 1-4 M Q ±AM NS , 
q = 3 ± Aq and for cl/Mbh = 0.9 ± 10% and 0.5 ± 10%. 
We see that the effect is small, but not negligible. 
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a)M NS =1.4 M sun ±AM NS 



1200 1400 
v GWtide ( Hz > 



q=5 q=3 




1200 1400 
v GWtide ( Hz > 



FIG. 6: (Color online) The same plot as in Fig.\^for a black 
hole with angular momentum a/ Mbh = 0.5 ± 0.05. The data 
for q = 10 and 15 are missing because no disruption occurs 
for the considered EOSs. 



An interesting feature common to the Figs. |31 E] and [5] 
is that, for the considered EOSs, all data are well fitted 
by a parabolic fit. 



It should be noted that the fit proposed by Vallisneri in 
ref. |H, i.e. 



Rns 



O.U5(i>M BH )- o n for vM BH < 0.045 
0.069(^M BH )-°- 95 for vM BH > 0.045 

(40) 

(where v = vawtide) for M NS = 1.4 M Q and M B h = 
(2.5—80) M Q , predicts values of Rns and VGWtide largely 
different from those we find, as shown in Fig. |8l In par- 
ticular we see that, for a given value of the NS mass, 
of the radius (i.e. of the stellar compactness) and of 
the mass ratio, the value of VGWtide evaluated in (l3j 
is systematically smaller than the value we find. This 
happens for a number of reasons. For instance, let us 
consider two binary systems having a black hole of the 



TABLE II: The values of q above which no disruption occurs, 
qmax, and the corresponding minimum value of vcwude, v-min, 
are tabulated for the binaries considered in Figs.\^and\Q i.e 
for a black hole with a/ Mbh = 0.5,0.9 and a neutron star 
with mass M = 1.2, 1.4 M©. 



a/ Mbh = 0.5 


EOS 


1.2 M 


1.4 M Q 




Qmax 


v min (Hz) 


Qmax 


Vmin (Hz) 


GNH3 


9.4 


780 


7.5 


840 


BGN1H1 


8.4 


878 


6.7 


942 


BPAL32 


8.2 


899 


6.6 


967 


BBS1 


7.6 


966 


6.1 


1032 


BPAL22 


7.6 


968 


5.8 


1084 


APR2 


6.9 


1068 


5.6 


1128 


BPAL12 


6.6 


1116 


4.6 


1375 


a/ Mbh = 0.9 


GNH3 


23.2 


559 


18.5 


596 


BGN1H1 


20.6 


628 


16.5 


667 


BPAL32 


20.1 


642 


15.9 


690 


BBS1 


18.8 


686 


15.1 


726 


BPAL22 


18.9 


689 


14.4 


759 


APR2 


17.0 


757 


13.8 


793 


BPAL12 


16.2 


790 


11.3 


960 



M NS= 1 4 M sun ± AM NS ■ Q=3±Aq 




1000 1200 1400 1600 1800 20I 

v GWtide ( Hz ) 

FIG. 7: (Color online) The plot of Rns versus VGwudc is 
done for the same NS mass, Mns = 1-4 Mq ± AMns, the 
same mass ratio q = 3 ± Aq, and two different values of the 
angular momentum, a — 0.9 ± 10% and a — 0.5 ± 10%. 



same mass and a neutron star with the same mass and 
compactness; in one case the NS structure is computed 
using the equations of Newtonian gravity and an n = 1 
polytropic EOS, in the other case using the TOV equa- 
tions and one of the EOS we consider in this paper. If we 
evaluate the value of r t ide for the first system using the 
affine approach in Newtonian gravity as in [l3j . and for 
the second system using our improved affine approach, 
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1200 1400 
v GWtide ( Hz > 



1600 1800 



FIG. 8: (Color online) The data of Fig.\3\are compared to the 
fit given in Eg. \40\ flj] - 



TABLE III: The parameters a, /3, 7 of our parabolic fits, a is 
in km, /3 in km/ Hz, 7 in km/Hz 2 . 









a/M B H 


= 0.5 








1.2 M Q 


1.4 M 


q 


a 


P 


7 


a 


P 


7 


3 


34.7 


-2.8-10" 2 


8.1-10" 6 


32.4 


-2.3-10" 2 


5.5-10" 6 


5 


38.5 


-3.6-10 -2 


11.9-10" 6 


35.3 


-2.9-10" 2 


7.5-10" 6 








a/MBH 


= 0.9 






3 


33.2 


-2.5-10" 2 


6.9-KT 6 


31.7 


-2. MO" 2 


5.010" 6 


5 


37.1 


-3.3-10" 2 


10.2-10" 6 


37.4 


-3.2-10" 2 


9.0-10" 6 


10 


50.0 


-6.7-10" 2 


27.7-10" 6 


56.3 


-8.0-10" 2 


34.5-10" 6 


15 


47.9 


-7.1-10" 2 


32.4-10" 6 


189.8 


-0.47 


31.4-10" 5 



We Will always find r tide Newtonian > r t ide improved- As a 

consequence the value of v GWtide evaluated by the New- 
tonian approach will be smaller than that evaluated with 
our improved approach. Furthermore, given the value of 
r tide , in [l3| fGwtide is calculated using the formula for 
a point mass moving on a circular orbit in Kerr space- 
time, while we compute this quantity using the 2.5 Post 
Newtonian equations describing the orbital evolution of 
a binary system. This introduces a further difference 
which makes our VQwude larger than that evaluated with 
the geodesic approximation, and the difference increases 
when we consider small values of the mass ratio q, as 
shown in Fig. [5J 

The parameters a, /3, 7 of our parabolic fits 



Rns 



■ ^VGWtide + lV Gwude 



(41) 



where Rns is expressed in km and vcwtide in Hz, are 
given in Table Hill 

We shall now discuss how the results of this paper 
could be used to estimate the radius of the star and to 
gain information on the equation of state, using gravita- 
tional wave detection. 




600 800 1000 1200 1400 1600 1800 20C 

v GWtide ( Hz ) 

FIG. 9: (Color online) Rns versus VGWtide is plotted for 
Mns = 1.4 M @ ±AM NS , q = 5±Aq anda/M BH = 0.9±0.09. 
The black continuous lines are the parabolic fits of the extreme 
points of the horizontal error bar for each EOS. The verti- 
cal dashed line corresponds to the supposed detected value of 
VGWtide- The two horizontal dashed lines identify the uncer- 
tainty on the evaluation of Rns- 



Let us assume that a gravitational wave signal has 
been detected and exhibits a frequency cutoff at, say, 
VGWtide = 1200 Hz. A suitable data analysis on the 
chirp part of the signal allows to determine the neutron 
star mass, the mass ratio, with the corresponding un- 
certainty given in Table HI and the black hole angular 
momentum with, say, an error of about 10%. Let us 
suppose that their values are Mns — 1-4 M Q ± AMns, 
q = 5 ± Aq and a/M BH = 0.9 ± 0.09, with AM ns and 
Aq given in Table [I] With this information, we compute 
the NS radius versus VQwude for the considered EOSs as 
explained in section III 21 and plot the results in Fig. |H1 
The two continuous black lines are the parabolic fits of 
the points which, for each EOS, are at the extreme of 
the horizontal error bar. If we draw a vertical line corre- 
sponding to the supposed detected cutoff frequency, i.e. 
VGWtide = 1200 Hz, we see that it intercepts the two lines 
in two points, which identify a region within which the 
stellar radius should fall. The data shown in Fig.[S]tell us 
that we would determine the NS radius with an error of 
2.5%. The error would be smaller if the detected value of 
VGWtide is smaller, and larger if it is larger. For instance 
for VGWtide = 1400 Hz it would be 3.3%, and 2.2% for 
VGWtide = 1100 Hz. 

Moreover, we would be able to exclude the EOSs which 
fall outside the region framed by the horizontal dashed 
lines, putting strict constraints on the equation of state 
inside the neutron star. 
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VII. CONCLUDING REMARKS 

In this paper we show that, by detecting a gravitational 
wave signal emitted by a BH-NS coalescing binary which 
exhibits a frequency cutoff due to the disruption of the 
star before the ISCO, we may be able to determine the 
radius of the star with quite a good accuracy, of the order 
of a few percent, and to put strict constraints on the 
equation of state of matter in the neutron star core. 

Our study does not intend to be exhaustive, since many 
more equations of state may be considered in the analy- 
sis; for instance we did not consider quark stars. There 
may exist other branches in the plots shown in Figs. |3][7] 
corresponding to more exotic EOSs. 

However, from our study it emerges that the quantity 
which mostly affects vcwtide is the stellar compactness; 
therefore, we expect that in general the cutoff frequen- 
cies lay on the parabolic fits which correspond to a given 
neutron star mass, mass ratio and black hole angular mo- 
mentum. 

A further point which should be discussed is the fol- 
lowing. In our analysis we have assumed that the value 
of VGWtide is known from the detection of a gravitational 
wave signal, but of course this quantity also is affected by 



uncertainties. For instance, we do not know how quickly 
the amplitude of the chirp goes to zero at tidal disruption, 
and therefore how sharp the step in the Fourier transform 
of the gravitational signal h{y) at vcwude is. To answer 
this question, an accurate modeling of the BH-NS coales- 
cence process is certainly needed. Moreover, assuming a 
given "slope" in h(v) at tidal disruption, how large would 
the experimental error in the determination of vcwude 
be? We plan to investigate this problem with a suitable 
data analysis study on the data of Virgo. 
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